Wall forces on a sphere in a rotating liquid-filled cylinder 
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We experimentally study the behavior of a particle slightly denser than the surrounding liquid in 
solid body rotating flow. Earlier work revealed that a heavy particle has an unstable equilibrium 
point in unbounded rotation flows. In the confinement of the rotational flow by a cylindrical wall a 
heavy sphere with density 1.05 g/cm 3 describes an orbital motion in our experiments. This is due 
to the effect of the wall near the sphere, i.e. a repulsive force (Fw)- We model Fw on the sphere as 
a function of the distance from the wall (L): Fw oc L~ A as proposed by Takemura and Magnaudet 
1 . Remarkably, the path from the model including Fw reproduce the experimentally measured 
trajectory. In addition during an orbital motion the particle does not spin around its axis, and we 
provide a possible explanation for this phenomenon. 



I. INTRODUCTION 

Forces acting on a sphere in various flows are of importance from a fundamental point of view and in applications 
(see the review article by Magnaudet and Eames [2 ). A natural approach for analyzing the forces in complex flows 
is to decouple the flow effect into those of typical flows, like a solid-body rotating flow. In recent years the forces on 
particles lighter than the surrounding fluid were studied in a solid-body rotation flow around a horizontal axis (i.e. 
gravity is perpendicular to the rotation axis) by Nierop et al. [3|, Bluemink et al. [HG3E], Rastello et al. [7, 8 . The 
advantage of this system is that buoyant particles reach an equilibrium point, from which drag and lift coefficient are 
well determined over a wide range of dimensionless parameters. Provided the equilibrium point is sufficiently far from 
the wall(s) bounding the rotating flow, this can be considered as unbounded. In this way drag and lift forces as well 
as particle spin have been measured at moderate Reynolds numbers and reported in the above mentioned papers. A 
particle heavier than surrounding liquid spirals outward in rotating flow, i.e., it has an unstable equilibrium point. 
Such orbits were calculated by Roberts et al. [9 for such small particles that Stokes flow can be assumed. They also 
took no lift forces into account. 

The effect of the wall is important but less intensively studied, especially at moderate-to-high Reynolds numbers. 
It has been investigated mainly in the situation of a sphere moving parallel to a fiat plate or rotating near a plane 
boundary. The influence of a vertical flat wall on a spherical rising bubble is studied experimentally by Takemura and 
Magnaudet [1 and numerically by Zeng et al. [10] for Reynolds number less than 100. Takemura and Magnaudet [1 
discussed two hydrodynamical mechanisms of wall interaction, one due to the vorticity generated at the bubble surface 
and the other to the irrotational dipole associated with the bubble motion. They proposed empirical correlations based 
on the strength of these two mechanisms, leading to practical expressions of the lift force as a function of Reynolds 
number and the distance between the wall and the particle. Liu and Prosperetti [11] studied the effect of the spin 
of the sphere in the presence of no-slip planar boundaries. They found that for small sphere- wall gap widths and 
Reynolds numbers ~ 1 viscous effect plays a role causing repulsive forces while larger Reynolds numbers ~ 100 cause 
a Bernoulli effect in the gap, which may turn the forces into attractive ones. 

When it comes to rotating flow, the wall effect is studied for very viscous flow. Ashmore et al. [12] investigated 
the interaction of the wall with a dense steel particle in the Stokes regime. Mullin et al. [13 found that with a very 
heavy sphere in a highly viscous rotating fluid three regimes can be distinguished as a function of the cylinder rotating 
speed, which we will discuss in § |IV| 

In this paper we experimentally investigate the trajectories of particles, slightly heavier than the rotating fluid 
in which the particles are introduced. In particular we are interested in the interaction with the cylinder wall for 
higher Reynolds number O(10 3 ). The translational and rotational motion of the spheres are observed with help of 
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FIG. 1: Forces acting on a particle: Fa the added mass force, Fd the drag force, Fl the lift force, Fg the gravity 
force or body force, Fj the force due to acceleration of the flow, and Fw the wall repulsive force. 



Particle Tracking Velocimetry (PTV) and the orientation detection software [14] , respectively. We use Particle Image 
Velocimetry (PIV) to record the velocity distribution in the undisturbed as well as in the disturbed fluid 

The outline of this paper is as follows. In § [IT] the equation of motion of a particle is introduced. The experimental 



setup is described in § III The results are shown in § |IV[ followed by conclusive remarks in §|V[ 



II. FORCES ACTING ON A PARTICLE 



The force balance on a sphere in the present condition is 

du 



P P V P 



dt 



Fa + Fd + Fl + Fg + -F/ + F w , 



(1) 



where p p is the density of the particle, V p the volume of the particle, u the particle velocity, Fa the added mass force, 
Fjj the drag force, Fl the lift force, Fg the gravity force or body force, Fj the force due to acceleration of the flow, 
and Fw the wall repulsive force. All forces in the rotating flow are shown in Fig. [I] The operators d/dt and D/Dt 
are material derivatives going with the particle and with the fluid, respectively. Owing to the fact that the present 
configuration is quite similar to that of light particles investigated by Bluemink et al. [6 , we can adopt expressions 
for Fa, Fd, Fl, Fg-, and Fj from Magnaudet and Eames [2 , Mazzitelli et al. [15 j. The full differential equation with 
all relevant variables is 



PpVp 



du 

~dt 



pVpC A 



f DU 
\~Dt 



du 

~dt 



1 



+pV p C L (U-u)x(VxU) + -pApC D \U-u\(U-u 



DU 

- [Pp ~ P) V p g+pV p — +F W , 



(2) 



where p is density of the fluid, Ca added mass coefficient (CU=l/2 for the sphere case), U the velocity of the fluid, 
Cl the lift coefficient normalized with the vorticity V x Z7, A p the cross sectional area of the sphere, and g the gravity 
acceleration. 

We have not found results on the magnitude of a wall induced force Fw in case of a curved wall, as is ours. Since 
the radius of the sphere is small with respect to the drum radius, we may confidently use results for a flat wall. Recent 
work, relevant for our case, is given by Takemura and Magnaudet [1 and Zeng et al. [ID] . From these work it is clear 
that there are two mechanisms which contribute to Fw- The first is the vorticity distribution in the wake behind 
the sphere. This diffuses outward, but this process is asymmetric due to the presence of the wall. It leads to a wall 
force away from the wall. On the other hand the accelerated flow through the gap between the sphere and the wall 
produces an attractive force. The first mechanism is dominant over a wide range of Reynolds numbers. Takemura 
and Magnaudet [I] did experiments at Reynolds number of order 10 2 and over a large range of the ratio between 
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distance from a wall and sphere radius, L/R. They summarized their results as follows. 

F w = C w A p pu 2 /2, (3a) 

where 

C w = C W0 (L*)a 2 (Re)(L/jRy^ (3b) 

with 

7 « 3.0, g(Re) « -2.0tanh (O.OliJe), a(i?e) « 1 + 0.6Re 1/2 - 0.55i?e ' 08 (3c) 

n (t*\-\ (9/8 + 5.78- 10- 6 )L* 4 - 58 /3 2 exp(-0.292L*) for < L* < 10 _ 
Cwo{L ) - | 8 94 ^2 L *-2.09 for io < L* < 300, [M) 

where /3 = 1 for a rigid sphere. The dimensionless separation L* = LU/v, the Reynolds number i?e = Ud p /v based 
on the particle diameter d p and the oncoming fluid velocity at the particle location U. 

Zeng et at. [10] performed numerical calculations for the situation of the experiments by Takemura and Magnaudet 
PP, and found good agreement with Eq. ([3^L-d) in the relevant parameter range. Although our Reynolds numbers are 
larger than those in the experiments of Takemura and Magnaudet [1 , we assume that the same trends with Reynolds 
number. 

We choose cylindrical coordinates (r, </>, z) with z along the rotation axis. Since during the experiments, to be 
described in the next section, the particle always stays in a plane perpendicular to the rotation axis, the momentary 
position is sufficiently described by r = re r + 0e^, where e r and are unit vectors in radial and azimuthal directions 
respectively. With a constant angular velocity cj, the flow velocity is 

U(r) =ujre 4> . (4) 

In terms of this coordinate system, Eq. ([2| becomes 

r = ^ - + CA ^ j {^Pp + Ca P) r ^ 2 + (P ~ Pp) 9 sin <t> ~ ^ r p C dP\J r 2 + (r (j> - u?j^ r + 2C L pur (u - <j?j 

- pruj 2 (C A + 1) + F w ) (5a) 
for the radial direction. For the azimuthal direction, we obtain 

4> = ( ^r^ ^ Cap) ) (~ 2 ^ pp + CaP " } ^ + ( p ~ p ^ 9 cos ^ + 8^ CdP \/ * 2 + ( r ~ u )) r { u ~ ^) + 2CL P uf ^j - 

(5b) 

The dots denote differentiation with respect to time. 



III. EXPERIMENTAL SET-UP 



Figure [2] shows the experimental setup. The cylindrical acrylic drum is 500 mm long and has a radius of 250 mm 
with 15 mm thick plastic walls and lids. Its axisymmetric axis is horizontal. Two steel rods with rubber coatings 
support it and one of them rotates driven by an AC servo motor. We vary the frequency of the drum fd in the range 
of 0-2 Hz. The drum is filled with de- ionized water with density 1.00 g/cm 3 . We use slightly heavier polysterene 
spheres with density 1.05 g/cm 3 . Their size ranges from 10 to 60 ± 0.01 mm. The adsorption of de- ionized water by 
polysterene is less than one procent by volume. 

Particle trajectories are recorded by using a high-speed camera (Kodak 2000, Redlake Co., USA / Photron SA.l, 
Photron Co., Japan) with a lens (Fujinon TV Zoom lens). Recording speeds are between 60 frames per second 
(fps) and 500 fps. In general, the measurement time is about 6 minutes which corresponds to more than 20 particle 
trajectory cycles. To ensure stable orbits, the drum is left rotating together with the inserted particle for more than 
45 minutes before starting each measurement with set drum frequency. 



FIG. 2: Experimental setup. 1. Laser, 2. Drum, 3. Drum controller, 4. Camera, 5. Camera controller 
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FIG. 3: Solid-Body rotational flow at a drum frequency of 0.60Hz: (a) A snapshot of the PIV tracer particles, (b) 
Flow field. The scratched region is highlighted in gray color, (c) Velocity vs. position on a horizontal line (red line 

in the left figure). 



The recorded images are analyzed with PTV software. The particle center detected by the use of a circular Hough 
transform method represents its position. The position of the wall of the drum is determined using particle reflection 
images. As shown in Fig. [3^a), the lid of the drum has scratched areas at the radial position of from 15.5 cm to 17 
cm. When the particle locates In the scratched area, we determine the particle position using linear interpolation. 

The flow in the drum is measured using PIV technique with a laser sheet (Lasiris Magnum 2, StockerYale on 
Coherent inc., Canada). The laser-beam illuminates a cross-sectional area of the drum. The flow filed at the drum 
frequency /^=0.60 Hz, without a sphere, is recorded at 1000 fps and an exposure time of 1 /is. Figure |3|b) shows the 



FIG. 4: Four regimes of a heavy particle trajectory: (1) Fixed-point regime, (2) cascading regime, (3) fixed solid 

body rotation regime, and (4) suspension regime. 



flow filed by averaging 250 frames. This shows a solid-body rotation. The velocity along the horizontal line (solid red 
line in Fig. [3Jd) is shown in Fig. [3^c) with the theoretical line of the solid body rotation. The scratched area of 15.5 
to 17 cm in radial position is indicated in gray color, where the mean velocity is not measurable. The experimental 
result agrees well with the predicted line within error bars, confirming that the flow is solid body rotation. 



IV. RESULTS AND DISCUSSIONS 



We describe what a particle, radius 7 mm, density 1.05 g/cm 3 does when immersed in the rotating flow and the 
drum frequency fd is increased from zero to 1.8 Hz. Much of its behavior is similar to that of heavy particles as 
described by Mullin et al. [13], mentioned in the introduction. For fd <0.07 Hz the particle rolls along the drum wall 
at a fixed position as regarded in the laboratory frame. This is what Mullin et al. [13] call the fixed point regime 
shown as trajectory (1) in Fig. |4j The other regimes mentioned by Mullin et al. [13] are found as well. For 0.12 
Hz< fd <1.2 Hz the particle touches the drum wall in a part of the orbit and falls down in the remainder. This is 
the cascading regime (trajectory (2) in Fig. H). For fd >1.2 Hz the particle sticks on the wall of the drum and rotates 
with it. This is the fixed solid body regime (trajectory (3) in Fig. [4|. 

However, we found that our particle, slightly heavier than the fluid, floats through the drum, without contacting 
the drum wall, at drum frequencies between 0.07 and 0.11 Hz. We call this the suspension regime shown as trajectory 
(4) in Fig. [4j Clearly this is possible by the presence of the wall. In the remaining part of this section we focus on 
the particle motion in this regime. 

The positions of the particle are plotted in Fig. [5] for several drum frequencies fd- The coordinates of the particle 
positions in Fig. [5] are made dimensionless with help of the particle radius. The trajectories show almost perfect 
circles. The minimum particles distance to the wall is ~ 1 for all cases. Their radius increases with increasing fd and 
their central positions (r e , e ), shown as cross symbols in Fig. |5j depart from the wall with fd- The coordinates r e 
and <j) e mark the equilibrium position. This position apparently is not stable in the sense that upon a displacement 
the particle returns to the original position. However, it does not spiral outward or inward as in the study of Roberts 
et al. [9 , but remains in orbit about (r e , <j) e ). The explanation for this lies in the fact that in the analysis of Roberts 
et al. [9 there is no lift force. At the equilibrium position (r e , (j) e ) in our case Cd and Cl can be calculated by taking 
r = r = <j> = (j) = in eq. [5] as: 



4 gd cos (fj e 

Cd= o 2~^~~ ( 6a ) 

3 riuj 1 



^ = -(1 + ^-^-^). (6b) 

From Eq. [5] it follows that in the r = r e , = situation the drag is vertical and balances effective gravity. Further 
the central force exerted by the pressure gradient in the rotating fluid is balanced by the lift force. When now the 
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FIG. 5: Particle position of a 7mm polystyrene particle in a drum for frequency 0.07 - 0.11 Hz. 




FIG. 6: Images of a slightly heavy particle with black and white texture in a drum 



particle undergoes a vertical displacement in r-direction, it meets with a drag which in excess of the buoyancy force 
and hence moves upward. Then, it acquires a 0, by which the lift is reduced and it is pushed inward by the pressure 
force. In the case of Roberts et al [9] there is no lift force and so the pressure gradient force is not balanced in r = r e , 
4> — 0. A picture of an orbit as measured by us is shown in Fig. [6] It is clear that between the points marked as 1, 2 
etc. the angular velocity varies slightly. It is also remarkable that, as the particle in the marked positions shows, the 
particle does not rotate about an axis through its own center. 

We are, of course, interested whether the Eq. ([5| are able to reproduce the observed orbits and the absence of 
spinning motion. The resulting range of the Reynolds number is 1400-1500. In this Re regime, the expression of 
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FIG. 7: The equilibrium radius vector r e , the momentary position of the sphere r, and the particle velocity u. The 
particle turns around r e with momentary angular velocity Q p . The dashed black line represents the particle circular 

orbit. The cross marker shows the equilibrium point. 



Eq. (I3p-d) can be reduced to g(Re) ~ —2.0, Cwo(L*) oc L 2 . Thus, the wall repulsive force may be written as 



F w w X L-\ 



(7) 



where x is a constant. We used in Eq. (5a) the expression ^ for the wall force Fw In addition we used x as a fitting 
parameter. We start with analyzing the trajectories. We denote the radius vector of the orbit center, measured from 
the center of the drum, with r e and of the momentary position of the sphere with r as shown in Fig. [7] . The particle 
turns around r e with instantaneous angular velocity fl p which varies along the orbit. So its velocity is 



u = Q p x(r - r e ). (8) 

This is the velocity as measured in the laboratory frame. The velocity relative to the surrounding fluid which 
rotates about the center of the drum with angular velocity u>, is 



u — Ljxr = (fl p — (jj)x(r — r e ) — u;xr e (9) 

The velocity —u)xr e is constant along the orbit and it directed along the negative y-axis in Fig. [5] We can make an 
estimate for this from the measurements of the orbits and the velocities. We take the case with the following data: 
R = 7 • 10 -3 m; fd = 0.11 Hz. The orbit is displayed in Fig.|8]and the velocity relative to the rotating frame in Fig. [9] 

We write Eq. ^ as 



u = —ujr e e y + Aftaee =< u > +u' , (10) 

where = \Q p — u>\ and a = \r — r e \. 

From Fig. [8] we obtain r e = 0.15 m and hence < u > should be 

< u >= 2nf d r e = 104 • 10" 3 m/8. (11) 

In Fig. § the maximum and minimum velocity are 111-10 3 m/s and 99-10 3 m/s, respectively and hence the 
average is 104 • 10 -3 m/s which agrees very well. From the difference of these values and considering only the ground 
frequency, that is 

u' = ACl max a cos (jut + 0), (12) 



where <p is a phase angle. We obtain u' = 6 • 10 -3 cos (uot + </>). From Eq. (10) it appears that the particle "sees" a 
constant velocity in the negative y-direction and a small fluctuation on top of that. For the calculation, with help of 
Eq. ([5| we therefore approximate Cd and Cl with the values based on < u > with 4> e =0. These are 

8 gR p,-p 
3 (wr e ) 2 p 



C L = \{1 + C A ). 



(13b) 
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FIG. 9: Time normalized by the drum time cycle vs. the magnitude of the particle velocity. 



With these values for Cd and Cl we calculated the particle trajectories for a number of frequencies in the suspended 
regime 0.07 Hz < fa < 0.11 Hz. We used the value of x to fit the calculated results to the measured ones. As 
comparison of the results in Fig. [6] with those in Fig. [8] shows, in this way an excellent agreement is found. The value 
of x is 5-10 -4 . 

Next we discuss the fact that the particle does not spin much around its axis. In order to make an estimate for the 
expected spin, which can only be induced by the fluctuation part of the velocity, we first look at the relaxation time. 
Bagchi and Balachandar [16] determined relaxation time for rotation of spheres in shear flow at Reynolds numbers 
up to 200 and found values about 4 times the quantity d 2 /v, where d is the particle diameter and v is the kinematic 
viscosity of the liquid. It is known that the relaxation time decreases with increasing Reynolds numbers. We have 
Reynolds numbers of order 10 3 and therefore estimate the relaxation time r as 



r = — - 200s 
v 



(14) 



9 



20 
15 
10 
5 



- absolute angle of particle 

- mean angle of particle 




14.2 



FIG. 10: The absolute angular rotation of the sphere in a solid body rotation flow with fd = 0.07 Hz for one period 

of the drum. 



This is much longer than the drum period, 

r » — . (15) 
UJ 

A model equation for the induced rotation is 

I0 + F0 = M, (16) 

where / is the moment of inertia of the sphere: / = 2mR 2 /5 with mass m = 47rp p R 3 /3. Then, dividing by /, F/I is 
1/t, and the equation becomes 

We have not found expressions for the moment in an oscillating flow. In general the moment exerted by the surrounding 
flow is expressed as 

M = \Cmp(u') 2 kR\ (18) 

Zeng et al. [10] report for a sphere in a shear flow and in the presence of a wall Cm values of 10 -2 -10 -3 . We take 
Cm= 0.01, insert this in Eq. ( [18] ) . Using also the above expression for / results in 

0+- ~l0- 2 e iujt . (19) 
r 

Taking = Be lUJt , we find for the real part of B 

10- 2 r 2 10- 2 

Re < B » - TT(^ ~ — ■ (20) 

With our uj = = 0.69 s _1 this amounts to an oscillation of 1.2 degrees. An oscillation of this order is observed 
indeed as Fig. [To] shows. 

As described above, the orbits from Eq. ^ with the wall repulsive force model can reproduce the experimental 
results. This suggests that even though the wall force model (Eq. [3| by Takemura and Magnaudet [1 was made for 
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FIG. 11: The result of the flow field around a 7 mm particle in a drum at /^=0.07 Hz using Particle Image 

Velocimetry. 



a rising sphere near a vertical planar wall, it seems to be valid for the case of solid-body rotation in the present Re 
regime. 



In addition, we performed PIV measurements. A snapshot of the flow field is shown in Fig. 11 'a). The inward flow 
is undisturbed and the wake of a particle is deflected inwards the cylinder. This is more clearly observed in Fig. 11 [b). 
The upstream flow at = —14.5° of the particle is undisturbed solid body rotation. The wake effect is observed in 
downstream of the particle at = 0.9°, = 6.3° and = 12.4°. 



V. CONCLUSIONS 



We study translational and rotational motion of a slightly heavy particle in a rotating drum filled with water. 
Remarkably in the regime of 0.07 Hz < fdrum < 0.11 Hz the particle is found to be suspended in the drum with an 
orbital motion, which was not observed by Ashmore et al. [12 j. We investigate the force balance by including the 
effect of the drum wall (Fw) modeled as Fw oc L~ A [TJ[T0]. The orbits of the particle from the force balance including 
Fw can reproduce the experimental trajectories. We also investigate the spin of a particle. The particle orientation 
changes less than 10 degrees over one cycle of the drum due to the small torque acting on the sphere. 
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